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Abstract. The electrical conductivity changes of the Si(lll)-(6x6)Au surface at early 
stage of Pb deposition was studied experimentally and theoretically as a function of 
coverage. Pb deposition onto a Si(lll)-(6x6)Au surface induce strong change of the 
conductance even at a low coverage (< 1 monolayer). The experimental results are 
analyzed using the theoretical model of two dimensional (2D) square atomic lattice. 

1. Introduction 

The electrical conduction through semiconducting surface layer is very sensitive 
to the presence of individual foreign atoms. These atoms act predominantly as 
dopants and can, depending on the semiconductor type, enhance or decrease the 
electrical conductivity of the depletion or accumulation layer pQ . Presence on the 
semiconductor surface of a monoatomic metallic layer allows to separate these 
quasi two-dimensional semiconducting effects from truly two-dimensional phe- 
nomena arising during interaction of individual atoms with the surface. In this 
case the conventional theoretical description which takes into account the statis- 
tics of current carriers in semiconductors or the well known Drude conductivity 
theory, is not sufficient. 

In this study we propose a theoretical model which allows to calculate the 
electrical conductance of truly two-dimensional lattice of atoms with single atoms 
deposited on it. The model simulates the initial stage of well known phenomenon 
of oscillating electrical resistivity observed during layer by layer growth of metallic 
layer [5). The model is illustrated by experimental data of the electrical conduc- 
tance of the Si(lll)-(6x6)Au system, measured during deposition of individual 
atoms of Pb on it. The morphology of the surface covered with Pb atoms is studied 
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with scanning tunnelling microscope (STM). 
2. Experimental setup and results 

The measurements were performed in a UHV chamber with a base pressure of 
less than lxlCT 10 mbar. The structure of the substrate and the deposition of Pb 
were monitored by reflection high energy electron diffraction (RHEED) system. 
An n-type Si(lll) wafer of around 25 flcm resistivity and 18x4x0.4 mm 3 size was 
used as substrate. Electrical conductivity was measured in situ by four-point probe 
method. An alternating current: 1 = 2 fj,A, 17 Hz was sent through the outer-most 
Ta clamps contacts, while the ac voltage was measured across the inner two W 
wires kept in clastic contact with the wafer. Before each measurement run, the 
surface was cleaned to obtain a clear Si(lll)-7x7 RHEED pattern, by few flash 
heating. In order to prepare the Si(lll)-(6x6)Au surface structure, 1.3 ML of Au 
were deposited on Si(lll)-7x7 superstructure. Annealing for 1 min at about 950 
K and slow cooling to room temperature (lOK/min.) resulted in the appearance 
of a sharp (6x6)Au superstructure RHEED pattern. During deposition of Pb the 
sample was kept at 78 K. The conductance could be measured simultaneously 
with deposition. The amount of deposited material in units of monolayer (ML = 
7.8xl0 14 atoms/cm 2 ) was monitored with a quartz crystal oscillator. 




Coverage (ML) 

Fig. 1: (a) The relative conductance changes recorded during Pb deposition onto a 
Si(lll)-(6x6)Au surface at 78 K. Go denotes the bare substrate conductance, (b) STM 
topographic image of Si(lll)-(6 x 6) substrate covered with 0.05 ML of Pb at 102 K. 
The smaller round bright species are the Pb single atoms. 

Fig. la shows the relative conductivity changes during Pb deposition at 78 
K onto a Si(lll)-(6x6)Au surface for coverage up to 0.2 ML. The decrease of 
the conductivity starts from the beginning of the deposition. The conductivity 
reaches a minimum at around 0.12 ML coverage and afterwards it increases. An 
abrupt increase of the conductivity at opening of the shutter is connected with an 
influence of radiation from evaporator. The samples studied with the STM were 
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prepared in similar condition in another UHV system equipped with Omicron 
VT STM. Pb films of various thickness have been deposited on such prepared 
substrates mounted in the cooled STM stage. During the deposition the substrate 
temperature was equal to about 100 K. The pressure of the chamber during the 
Pb deposition was below 7 x 10 11 mbar. Fig. lb shows STM scan taken after 
deposition of about 0.05 ML of Pb. Previous STM studies revealed [H] that at 
this coverage the visible small features are predominantly single Pb atoms. A few 
larger and brighter species are probably composed of more than single atom. The 
density of these atom clusters increased with Pb coverage and finally, for coverage 
corresponding to 1 ML of Pb(lll), the whole surface was covered with epitaxial 
monoatomic layer of Pb. 



3. Theoretical description and results 

In this section we study theoretically the conductance through a model of two- 
dimensional (2D) square atomic lattice. The model is schematically shown in 
Fig-EK- The central region, between the left (L) and right (R) electron reservoirs, 
consists of an array oiNxM atom sites. Additional atoms (adatoms) are added on 
the top of 2D lattice and cause the disturbance of the electron transport through 
the system. A group of coupled adatoms can be treat as a cluster of atoms. The 2D 
lattice of atoms is considered as the effective substrate which corresponds to, used 
in the experiment, Si(lll)-(6x6)Au substrate. The additional atoms correspond to 
deposited Pb atoms and the left and right leads represent the voltage electrodes. 
We assume here, that each adatom is coupled only with the single substrate atom. 
Moreover, 3D island growth is not considered. 

The total tight-binding Hamiltonian for the system can be decomposed as 
H = H + Hi + Vi where H describes electron states in the left and right 
electrodes, in the substrate and also interactions between substrate atoms. The 
adatom single electron states are described by the Hamiltonian H\. V\ concerns 
adatom-adatom and substrate atom-adatom interactions. For simplicity, a single 
orbital per atomic site is considered. Only nearest-neighbor interactions are taken 
into account. We consider the total Hamiltonian H in standard second-quantized 
notation in the following form: 
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N I\I N M-l 

Vi = Vdc ij ai J A v v dm cf j c i:j+1 A ij A ij+1 

i=l j — 1 i—1 j — 1 

N-l M 

+ ^ ^ VdncjjCi+ijAjj A i+ ij + h.c, (2) 

and Hi = J2iLi J^jLi £ tj c tj c ij A ij ■ Here the operators o-k a (a^ a ), are the annihi- 
lation (creation) operators for the reservoirs conduction electron with the wave 
vector k, (a = L,R). aij(o^) represents the annihilation (creation) operators for 
the electrons localized on the substrate atom ij (i = 1, N, j = 1, ...,M) and 
Cij (cfj ) on adatoms above ij substrate atom, and efj are the single electron en- 
ergy levels of the substrate atoms and adatoms, respectively. The elements V^ L ^ 
are the hybridization matrix elements responsible for electron transport between 
corresponding substrate atoms and electrodes; V n i m ) between the lattice atoms; 
Vdn(dm) between adatoms and Vd - between adatom and the 2D lattice atom. The 
function Ay in the Hamiltonian is equal to 1 if there is an additional atom at 
site ij. Otherwise it is equal to zero. This function allows to control number of 
adatoms included into calculations. Note, that the first term in Eq. 2 describes 
the interaction between the adatom and the substrate atom at site ij. 



(a) (b) 




Fig. 2: (a) The top view of the system of two-dimensional square lattice of atoms (N x M) 
connected to the left (L) and right (R) leads. The additional atoms are deposited on the 
top of this lattice; (b) Side view of 2-atom chain with two adatoms. V n , V m , Vl, Vr, Vd 
and Vdm represent the hybridization matrix elements responsible for electron transport 
between corresponding atoms. 

Electron and transport properties are analyzed within the framework of the 
Green's function technique. In the energetic representation Green's function sat- 
isfies the following equation of motion 0] : 

EG r ^ 2 (E) = {[a Pl -a+ 2 ] + ) + ([{a 01 ,H];a+ 2 }) E (3) 

Here E is the energy and /Jj represents the corresponding atomic site. The knowl- 
edge of the appropriate Green's function allows to obtain the electron local density 
of states (LDOS) of the substrate atoms or adatoms, i.e. for the site /3 one can 
write: LDOSp — ~^ImG r p p(E). The zero-temperature linear conductance, ob- 
tained for zero or very small voltage, is given in terms of the Green function by the 
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following expression: G{E) = ^-T(E), where T{E) = ^Tr{r L G r r R (G r )*} is 
the transmittance of the system. The matrix Yl is expressed as follows: (T^ij^i — 
T L S a 5 kl , where T L = 27r£g L \Vj: L \ 2 S{e - £fc L ) = 2ttVI/D. Here D is the effective 
bandwidth of the left lead density of states. The similar expression we can write 



for T R matrix, i.e. (T R ) ijtk i 
of the system becomes: 



F R 5iM3kM ■ Using above relations the conductance 
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G(E) = ^—T L T R 



N M 



G li,Mj 
i=l j=l 



(E) 



(4) 



Note, that the conductance is expressed only by the Green's functions Gl iM j(E) 
between atoms in the first and the last column, cf. Fig. (adatoms are not 
connected to the leads). The appropriate Green functions are obtained using Eq. 
3 and the Hamiltonian. The general matrix equation for G£ ■ kl (E) can be written 
in the form: A ■ G r = I, where / is the unit matrix. For the substrate atoms the 
matrix A is expressed as follows: 



A-ij^ki{E) — (E — Sij)Sik5ji 

—V m (Sik5j + u + 5ik5j-u) 



pL pfi 



(5) 



V n (Si +lk Sji + 5i-ik5ji) - V d Ai 



and similar for the adatoms. The retarded Green functions are obtained by finding 
the inverse of the matrix A, i.e. kl {E) = Aj^^E). 

Next, we apply the above formalism to the case of two-atom lattice connected 
to the left and right leads and two adatoms as shown Fig. \^p. This simple model 
allows to study the influence of the couplings between adatoms and also between 
these atoms and the substrate atoms on the LDOS and conductance through the 
system. We treat the Fermi energy of the system as the reference energy point 
i.e. Ep — (the linear conductance is obtained for the left and right chemical 
potentials equal to the Fermi level, fiL — M-R — Ep), and take into consideration 
the case of T L = T R = T. Moreover, the same electron energy of all substrate 
atoms and adatoms are assumed i.e. = eo ; ^tj = £ d- This assumption is quite 
reasonable as only linear regime is considered. All energies are given in T units 
and the conductance in 2e 2 /h units. 

The transport properties of the system are dependent on the density of states 
of this system. The conductance is proportional to the DOS of the substrate at 
the Fermi level 0] . For the considered system the matrix A can be written in the 
form: 



A = 





- £ + iT/2 


-Vm 


-Vd 


\ 




-v m 


E-e + iT/2 





-v d 




-v d 





E-e d 


— Vdm 
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-v d 


— Vdm 


E-e d / 



(6) 



By calculating the inverse of the matrix A we find the appropriate Green's function 
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needed to obtaining the LDOS. For the sake of symmetry of the system the 
LDOS of the substrate atoms and adatoms are equal, i.e. LDOSi — LDOS2 and 
LDOS3 = LDOSi - the indexes 1 and 2 (3 and 4) correspond to the substrate 
atoms (adatoms). The influence of the coupling between adatoms, Vd m , on the 




Fig. 3: The LDOS of the substrate atoms for the system schematically shown in Fig. 
^Jp. The couplings between atoms are: Vd = 2, V m = 2 and Vdm = 0, 1.5, 3 - solid, dotted 
and doted-thick lines, respectively. The dashed line corresponds to the case: V m = 2, 
Vd — (no interaction between adatoms and the substrate atoms), T = 1. 

LDOS of the substrate atoms is presented in Fig. [31 The dashed line corresponds 
to the case Vd = - no interaction between adatoms and the substrate atoms. 
The others curves correspond to the coupling between adatoms: Vdm — 0, 1.5,3 
- solid, dotted, dotted-thick lines and Vd — V m = 2. The electron energy levels 
are £0 = and £d — 0.5 (we assume that the adatom level, in comparison to the 
substrate-atom level, possesses the higher energy because the ionization energy 
of Au is greater then for Pb atom). The dashed curve possesses two maxima for 
E ~ ±V m which is known behaviour of atomic wires e.g. [5]- The additional atoms 
blockade the LDOS around the adatom energy level E = Ed, cf. the solid line. For 
increasing value of the coupling Vd the LDOS dip broadens (not shown here). 
However, for nonzero coupling strength between adatoms the LDOS increases in 
the regime of energies around the Fermi energy (so about Ed and £q) - cf. the 
dotted lines in Fig. El It causes that the conductance of the system increases with 
increasing V dm - 

The main conclusions from this simple model are: (i) non-coupled adatoms 
on the 2D lattice cause decreasing of the LDOS of the substrate (and also the 
conductance) obtained for the Fermi energy while (ii) the adatoms which are 
coupled with each other cause increasing of the LDOS at the Fermi level (so the 
conductance through the system is greater then for non-coupled atoms in the 
cluster). 

Now we take into consideration the system of the two-dimensional substrate 
of N x M atoms and adatoms deposited on it. A series of calculations were 
performed for the system with increasing number of adatoms. These adatoms are 
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Fig. 4: The relative conductance for the lattice of atoms: N = M = 12 (left panel), 
14 (middle panel), 20 (right panel) versus the coverage of the substrate. The other 
parameters are: V d = 2, V n = V m = 3, V dn = V dm = 5 and e = 0, e d = 0.5. 

shows the relative conductance G/Go (Go is the conductance of the substrate 
without adatoms) for the lattice of atoms: N = M = 12 (left panel), 14 (middle 
panel), 20 (right panel) versus the coverage of the substrate. The lines are plotted 
for better visualization and only the coverage of order 20% is considered as the 
time computation increases very strongly with increasing number of atoms in 
the system. For greater number of atoms in the 2D lattice the greater number 
of points is observed (each cross corresponds to one adatom deposited on the 
substrate). At the beginning the conductance decreases with increasing coverage 
of the substrate. This is due to the single atoms deposited on the lattice (cf. 
also Fig.|2J). However, for the substrate coverage about 0.12 ML, the conductance 
possesses the local minimum. In that point the first cluster of adatoms appears 
and for greater coverage more clusters are observed. These clusters are responsible 
for increasing the conductance although there is no percolation between the left 
and right electrodes. This effect is also visible in Fig. |3| for the one-dimensional 
case. Note, that the results presented in Fig. 0] could be somewhat different for 
various distribution of the addatoms on the substrate. But it was checked that 
the general results remain unchanged. 

To conclude, the theoretical results show that the conductance of 2D atomic 
lattice decreases with increasing coverage. For greater coverage (> 0.1ML), due to 
clusters of adatoms, the local minimum in the conductance appears. These results 
are in good agreement with the experiment described in the second section and 
also with the STM measurements. Due to the simple theoretical model the results 
presented in this section should be considered only as the qualitative conclusions. 
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